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ABSTRACT 

Red quasars are a population, characterized by significant extinction in UV, which could be explained by ab- 
sorption of dusty gas on a scale of a few kpc. We show that the enhanced radiation-pressure drives the dusty gas 
to supersonically expand and produces shocks. The shocks energize electrons to be relativistic via the first Fermi 
acceleration. As a balance of shock acceleration and synchrotron emission and inverse Compton scattering, the 
maximum Lorentz factor of the electrons reaches as ~ 10 6 . The shocked interstellar medium appears as extended 
multiwavelength fuzz, in which synchrotron emission from the electrons peaks at near infrared or UV bands and 
inverse Compton scattering around LOGeV-O.lTeV Future multiwavelength images of the fuzz would provide 
new clues to study the details of radiative feedback if red quasars could be a certain phase in evolutionary chains 
of galaxies. 

Subject headings: black hole physics — galaxies: evolution — quasars: general 



1. INTRODUCTION 

Radiative feedback from quasars to surroundings of their 
hosts has been realized for many years. Several effects have 
been suggested, including the Compton heating (Begelman 
1985; Chang et al. 1987; Wang et al. 2005), acoustic wave heat- 
ing (Chelouche 2007), proximity effects (Bajtlik et al. 1988; 
Goncalves et al. 2008; Wild et al. 2008; Prochaska & Hen- 
nawi 2008). The strong feedback due to radiation may quench 
the inflows fueling the central black hole (e.g. Binney & Ta- 
bor 1995) and star formation so as to arise the relation between 
black hole masses and dispersion velocities (Silk & Rees 1998; 
King 2003; Ciotti & Ostriker 2007; Fabian et al. 2006). How- 
ever, the direct large scale evidence for feedback is scarce, and 
even absent in most cases. 

Traditionally quasars are selected based on their bluer color, 
however, the existence of significant population of red quasars 
has been the subject of debate for a number of years. The 
percentage of red quasars is only ~ 6% or so from the large 
database of Sloan Digital Sky Survey (SDSS) (Richards et al. 
2003), - 20-30% in the FIRST-2MASS red quasar survey 
(Glikman et al. 2007), and ~ 20-30% in Spitzer survey (Brown 
et al. 2006; Lacy et al. 2007). The fraction of red quasars re- 
mains a big debate, however this might be due to heterogeneous 
definition or selection effects of surveys. There is increasing 
evidence for dust extinction as one of the major explanations 
for red quasars from SDSS (Richards et al. 2003), GALEX data 
(Trammell et al. 2007), Spitzer survey (Brown et al. 2006; 
Lacy et al. 2007) and HST images (Urrutia et al. 2008) as well 
as from numerical simulations (Di Matteo et al. 2005; Hop- 
kins et al. 2007). The debris, ejected from mergers or starburst 
triggering the central black holes, surrounds the nucleus prior 
to the galactic winds (Sanders et al. 1989; Yun et al. 2004; 
Fabian et al. 2006 and Hopkins et al. 2007), reddening quasars. 
The well-known roles of dust particles in outflows, as results 
of strong coupling between the charged grains with gas (e.g. 
Draine & Salpeter 1979), have been realized in the dynamics 
of dusty gas (Chang et al. 1987; Murray et al. 2005; Fabian et 
al. 2006). It is then expected that the swept medium by quasar 
radiation results in extinction of quasars. 



The present Letter makes an attempt to connect the action of 
the radiative feedback to the phenomena of red quasars and pre- 
dicts some observational features. We argue that radiation from 
quasar is partly converted to kinetic luminosity of supersonic 
outflowing gas and then partly tansformed to the kinetic energy 
of accelerated electrons via Fermi mechanism, producing ex- 
tended fuzz with non-thermal emission from radio to 7-rays 
around red quasars. 

2. DYNAMICS OF DUSTY GAS 

Dusty winds in homogeneous medium have been discussed 
in the literatures (Chang et al. 1988; Murray et al. 2005; Fabian 
et al. 2006). On the other hand, the gas in dark matter halos is 
driven into two phases in temperatures by the thermal instabil- 
ity (Fall & Rees 1985), especially in the presence of radiation 
heating of quasars. Thermal states of the medium are not the 
main goals of this Letter, but we will discuss two cases of the 
medium: 1) homogeneous interstellar medium (ISM) and 2) 
clumpy ISM, respectively. Observational appearance of the ex- 
tended fuzz is quite similar for two cases unless the high spatial 
resolution of radio telescopes could resolve the cold clouds. 

2.1. Homogeneous ISM 

For an isothermal sphere of gas with density profile as p oc 
R~ 2 , Murray et al. (2005) show the velocity of the momentum- 
driven winds at R for the optically thin limit 
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from eq. (26) in Murray et al. (2005), where R is the ra- 
dius of the expanding winds, R exp = o^a-r /2TiGm p « 123 cr9( )0 pc 
is the initial radius of the winds (Fabian et al. 2006), a is 
the dispersion velocity (a = 2OOCT200 km s" 1 ), Lndd = 2.07 x 
10 46 <T2oo er S s S_I i s tne Eddington luminosity, and <tj is the 
Thompson cross section. 

Here we only consider graphite grains here and assume they 
are sphere with a radius a. We assume the ISM has the same 
gas-to-dust ratio with the Galaxy. Since the charged parti- 
cles of dust are strongly coupled with the ionized atoms via 
Coulomb interaction (Draine & Salpeter 1979), the dusty gas 
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feels a stronger radiation pressure, amplified by a factor of 

A= — = 4.7 x 10 2 a 2 5 , (2) 

where cr d = 3.1 x 10~ 22 a 2 5 cm" 2 is the mean cross section per 
hydrogen atom for UV photons in the presence of dust (Martin 
& Ferland 1980; Chang et al. 1988; Fabian et al. 2006) and 
a_5 = a/10 -5 cm. The factor A has been discussed for various 
cases in detail by Fabian et al. (2006). For dusty gas, the Ed- 
dington limit in eq. (1) should be replaced by L' Edd = Lem/A. 
For a quasar with a bolometric luminosity LboI, the UV lumi- 
nosity is approximate to Luv = C^Boi = C^Edd, where £ w 0.1 
(Marconi et al. 2004) and A is the Eddington ratio. Inserting 
Lg dd into eq. (1), we have the maximum velocity of the expan- 
sion 

Vmax = 2 ( r v /^AC-ln(^AO-l« 13-4(7.4^5 (A£ .i) 1/2 , (3) 
at the radius 

fl w = A\£R exp » 6.0 ^50oA£ .i 4)0 kpc, (4) 
and the Mach number is 

M m - dX = V m Jc s = 93 .4 X£ (A£ .i) 1/2 T~ l/2 , (5) 

where ^ 5 oo = .4/500, £ .i = Z/0.1, c s = 2.87 x 10 6 r 5 1/2 cm s" 1 
is the sound speed and T$ = T/10 5 K. Here we assume warm 
winds with a temperature of 10 5 K. The supersonic winds form 
shocks inevitably and the strongest shocks are at R v w 6.0kpc. 
It should be noted that the velocity from eq. (1) is insensitive 
to the radius after R v , implying that strong shocks still ex- 
tend further. Chang et al. (1987) showed that shocks appear in 
both one- and two-dimension flow by detailed dynamical cal- 
culations, which agrees with the present simple estimation. As 
a result of so strong shocks, electrons will be accelerated by the 
first order Fermi process (Blandford & Eichler 1987). 

Given the extinction coefficient Ay, the fraction of radi- 
ation absorbed by the dusty gas is a = 1 - 10~ 04Av . We 
assume that the absorbed energies are converted into the 
kinetic luminosity of the winds L^n = ^C^Edd = 1.8 x 
10 44 77o.5ao.nA£o.i02 00 ergs s~\ where ao.17 = a/0.17 for Ay = 
0.2 and 770.5 = r//0.5 is the processing efficiency. The total en- 
ergy channeled into relativistic electrons is given by L n on-th = 
/ a Lkin, where /„ is the acceleration efficiency converting the ki- 
netic energy into the electrons. Though the detailed process of 
the acceleration is still insufficiently understood, it is generally 
taken as f a = 0.05 (Blandford & Eichler 1987). We then have 
the non-thermal luminosity from the electrons 

inon-th = 0.9 x 10 43 / 05 770.5 ao.nACo.icrjoo er S s s_1 ' ( 6 ) 
where /0.05 = / a /0.05. The spectrum emitted from the electrons 
will be given in §3. 

2.2. Clumpy ISM 

The physics of the infalling gas has been extensively dis- 
cussed stemming from Fall & Rees (1985). The basic con- 
straints on the medium are from the gravitational collapse and 
thermal conduction between the hot and cold phases (Fall & 
Rees 1985; Mo & Miralda-Escude 1996). The typical tem- 
perature and density of the hot gas are of 7), = 10 6 K and 
iif, = 10~ 2 cm" 3 in elliptical galaxies, respectively (Mathews & 
Brighenti 2003). 

The lower limit of cold clouds is given by M con d = 5.0 x 
lO 4 /?g kpc M due to thermal conduction, otherwise they will be 
evaporated (Fall & Rees 1985), whereas the clouds will collapse 



into smaller ones if they are too massive (e.g ~ 10 6 M Q ). We as- 
sume the clouds are sphere with radius R c and mass M c = qM con a 
(q > 1). A pressure balance n c T c = nhTh holds between the 
clouds and their surroundings. The geometric covering factor 
can be defined as 

where Aft is the solid angle of an individual cloud and Af c is 
the total number of the cold clouds. This covering factor is the 
same with the probability that a quasar is a red one. Here we as- 
sume that all the clouds are identical with same size, shape and 
physical conditions. For simplicity, we assume that the clouds 
are sphere. If the radius of the clouds is R c and R c <C R, we 
have Afi /4ir = R 2 /2R 2 , where R is the distance of the clouds to 
the galactic center. The number of the clouds is then estimated 
from eq. (7) if the geometric covering factor is estimated from 
the percentage of red quasars. The individual cloud should not 
be too dense, otherwise the red quasars become type 2s. This 
constraint can be set up by the column density tV h through the 
extinction coefficient Ay. We then have the density and temper- 
ature 

n c = 8 A q~ 1/2 N% 2 cm' 3 , (8) 
T c = l.2xl0 3 q^ 2 n h ,- 2 T hi6 K, (9) 

the radius R c = 38.7 q l / 2 N 2 ^ 2 pc and the cloud number is 
AC = 1.3 x 10 3 C R o.iq- l N 2l Rl k]K , where C R0 .i = C R /0.1, T hfi = 
7/,/10 6 K and nu-i = n/,/10 -2 cm" 3 . Though we do not know 
the details of the parameter q, we regard eq. (8-9) as the typical 
values for n c and T c since they are not very sensitive to q. The 
Mach number of the dusty winds in the clouds is 

M = 9.3 x 10 2 ^ 2 (ACo.i) 1/2 r 3 - 1/2 , (10) 
where T3 = r e /10 3 K. The Mach number is stronger than the ho- 
mogeneous ISM since the temperature is lower in the clouds. 
The shocked clouds produce relativistic electrons. 

When a cloud is just on the line of an observer's sight, the 
quasar appears as one red quasar. The transparent luminosity 
through the clouds is given by L\ = eAf c (Af2/47r) L m , where 
e = I0~ a44v and L; n is the intrinsic luminosity of the quasar, 
while the unabsorbed luminosity is given by Li = (1 -CR)Lj n . 
The global energy budget is then given by L obs = L\ +L 2 , we 
have L obs /L in = eC R + ( 1 - C R ) and 1 - L obs /L in = ( 1 - e)C R means 
the fraction of the intrinsic luminosity converted into the ki- 
netic energy of the winds. The kinetic luminosity of the winds 
is then given by Lkin = (1 - e)CR?7^ALEdd- The total non-thermal 
luminosity from the relativistic electrons is given by 

£non-th = / fl (l-e)C R 77£AL Edd . (11) 

Eq.(l) is still valid for the clumpy ISM (Murray et al. 2005), 
but the shocks are produced mainly on the cold clouds. 

3. ELECTRON ACCELERATION AND EMISSION 

Some electrons will be energized to be relativistic through 
Fermi shock acceleration. The maximum energy of the rel- 
ativistic electrons is determined by the energy gain and loss. 
Here we focus on the non-thermal radiation from the relativis- 
tic electron, rather than the detail of acceleration process. In 
this Letter, we assume the synchrotron radiation and inverse 
Compton scatter are two main non-thermal radiation processes. 
Since the magnetic field is poorly understood in such a con- 
text, we assume an equipartition between the tangled magnetic 
field and the gas, then the energy density of the magnetic field 

Ub = 1.38 x 10- 12 (« e r c ) 4 ergs cm" 3 and B = 5.9 (n c T c )\ /2 fiG, 
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where (n c T c )4 = « c r c /10 4 K cm" 3 . There are three sources of 
the seed photons of Compton scattering: 1) from the nuclear 
of quasars; 2) the reprocessed photons by the clouds; 3) syn- 
chrotron photons and 4) from the host galaxies. However, light 
from hosts can be neglected for quasars bright enough. The 
second could be neglected since they are only ~ 10% of the 
first source (if the optical depth r « A v for A v = 0.1, but may 
be important for the very red quasars or redden quasars). The 
energy loss is mainly due to synchrotron radiation and inverse 
Compton scattering of the seed photons from quasar's nuclear. 

The energy density of quasar's radiation is given by U qso = 
7.7 x 10" n L46/?6 kpc ergs cm -3 at R from the center, where 
R 6kpc = fl/6kpc and L 46 = L/10 46 ergs s" 1 . We find that U qso > 
Ub, indicating that the energy loss is mainly due to external 
inverse Compton scattering provided by the red quasars. The 
synchrotron self-Compton scattering can be neglected in such 
a context. For simplicity, we assume the quasar spectrum is 
characterized by a black body with a temperature of 10 4 K from 
the standard accretion disk (Inoue & Takahara 1996). The 
timescale of energy loss due to the inverse Compton scatter- 
ing IS fi oss = 3.0 x loV'tr 1 sec, where U w U qso , whereas 
the acceleration time scale is f a cc = R L c/v 2 h (Blandford & Eich- 
ler 1987), where Rl is the Lamore radius, and v s h ~ Vmax ~ 
10 8 cm s" 1 is the shock velocity. The maximum Lorentz factor 
of the electrons is given by fi oss = t acc , 

7max = 7.7 x 10 6 v shi sU:ll 2 B!i 2 , (12) 
where i/_io = J/ q s O /10~ 10 ergs cm -3 , v s h,8 = v s h/10 8 cms _1 and 
B-(, = Z?/10~ 6 G. The synchrotron radiation gets a peak at 

^s y „=4.5xl0 13 B_ 67 L,7Hz, (13) 

where 7 max j = 7max/10 7 . The peak frequency due to inverse 
Compton scattering is 

He w 7 2 ^disk = 3.6 x 10 27 ( 7 /7max) 2 T 4 Hz, (14) 
where = 3kT^/h, h is the Planck constant and 7dj s k = 
10 4 T4K is the maximum temperature in Shakura-Sunyaev disk 
for quasars. We assume a simple power law of the electrons as 
N{"f) = Noj~" in a range of 7 m ; n < 7 < 7 max via the first order 
Fermi acceleration, where j m \ n is the minimum Lorentz factor. 
We treat 7 m i n as a free parameter in this paper. iVo is determined 
by eqs. (6) or (11). For a red quasar containing a black hole 
with M, = 10 8 M Q , we plot the multiwavelength continuum of 
the fuzz for three different n in Fig 1 . It should be noted that 
the difference between v\c and the spectrum in Fig. 1 is due to 
the inclusion of the Klein-Nishina effects in the later. 

Some uncertainties of the proposed model should be noted, 
which are most likely because of the dependence of the factor 
A on the thermal state of the medium as shown in Fabian et al. 
(2006). Here we use the maximum A. The 7 max , depending on 
the shock velocity, could be lowered when A decreases, leading 
to the synchrotron peak shifts to long wavelength as well as the 
inverse Compton peak. We use a magnetic field in an equiparti- 
tion with thermal energy of medium, but it can not be deviated 
from the equipartition too much. 

4. DETECTING THE FUZZ 

The isothermal sphere under the galactic potential gives 
the density p{R) = f g a 2 /2irGR 2 , we have the column den- 
sity of the gas exterior to R exp , Nh = f g a 2 /2TrGm p R = 3.0 x 

10 21 /g0.icr2oo^6k P c cm ~ 2 ' where fg = -!/g,o.i is the g as frac- 
tion (Fabian et al. 2006). This leads to an extinc- 
tion coefficient A v = 0.2QJq A^2i, appearing as a red quasar, 



where N 2 \ = 7V H /10 21 cm" 2 and £500 = Q/500 is the gas-to- 
dust ratio (Gorenstein 1975). We note that the sweeping 
timescale of the radiation-driven winds is r w = R/V(R) = 2.2 x 

10 6 /? 6k pc (V/V max )" yr. The lifetime of quasars T qso is about 
the Salpeter time T Sa i p = 4.5 x 10 7 r7_o.i yr (77 = ?7/0.1 is the ra- 
diative efficiency), we have the percentage of red quasars at a 
level of r w /T qso ~ 5.0%. This roughly agrees with the results 
from SDSS (Richards et al. 2003). 
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FIG. 1 . — The multiwavelength continuum from radio to 7— ray bands. We 
take L non _ th = 10 43 ergs s" 1 , 7max = 7.7 X 10 6 , 7min = 900, U = 1QT 10 ergs cnT 3 
and B = 6.O/1G. The solid, dotted and dashed lines correspond to the indexes 
of n = 2.5, 2.7 and 2.3, respectively. The Klein-Nishina effects are included in 
the inverse Compton scattering. 



Table 1. Detection of the multiwavelength fuzz 
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note 
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ALMA 


47.4 


16.8 


yes 


0.09 


0.3 


yes 


MIPS 


3.0 


1.0 


yes 






no 


JWST 


0.4 


0.14 


yes 


0.7 


2.6XIO- 3 


yes 


GLAST 


150 




yes 


0.3 




yes 


HESS (I) 


0.03 




marginal 






no 


VERITAS 


0.03 




marginal 






no 


HESS (II) 


0.03 




yes 






no 



Note: 8 is the angular radius (in arcsec) at z. Col. (1) is the name of instru- 
ment; Col. (2-4) flux (in mJy), brightness (/dy arcsec -2 ) for ALMA, MIPS 
and JWST and note, respectively, but in 10~ 9 phs s~' cirT 2 for 7-ray detectors; 
Col. (5-7) are the same with (2-4), but in /dy, /dy arcsec -2 for z = 0.2. "yes" 
stands for detection, "no" for non-detection and "marginal" marginal detection. 
We assume the cosmological constants H Q = 71 Mpc~' km s~', S7m = 0.27 and 
Q A = 0.73. 

ALMA: Atacama Large Millimetre/Submillimetre Array, spatial resolution 
A9 = 0.005" at v = 900GHz. 

MIPS: Multiband Imaging Photometer for Spitzer, A6* = 5", working at 24, 70 
and 160/im. Here is for 24/i. 

JWST: J. Webb Space Telescope. Here is for NIR camera at 5/^m. 

GLAST: Gamma-ray Large Area Space Telescope (> lOOMeV), sensitivity 

~ 10 -9 phs s~' cirT 2 . 

HESS: High Energetic Spectroscopic System (> lOOGeV), sensitivity ~ 10~ 10 
(Phase I) and ~ 10" 1 1 phs s~' cnT 2 (Phase II). 

VERITAS: Very Energetic Radiation Imaging Telescope Array System 

The targets for searching the fuzz can be found from those 
low redshift red quasars from FIRST-2MASS, Spitzer and 
Sloan samples. Table 1 gives the future detections of the fuzz 
through multiwavelength telescopes. The predicted fuzz has a 
level of synchrotron power of 10 40 ~ 41 ergs s" 1 , which is easily 
to be detected by radio telescopes. IR band telescopes are obvi- 
ously capable of imaging the fuzz, in particular, the ALMA and 
JWST are able to detect the fuzz at z > 0.2 in future. The fuzz 
could be marginally detected by the existing HESS and VERI- 
TAS, but it is worth making an attempt. GLAST is anticipated to 
explore fuzz more feasibly. We have to note that TeV emission 
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from fuzz could be diluted via 7-7 interaction by the cosmic 
IR back ground photons beyond redshift z = 0.2. So it is not 
expected for 7-ray telescopes to detect the TeV emission from 
fuzz for red quasars with z > 0.2. If the episodic activity of 
supermassive black holes is really caused by a series of merg- 
ers of galaxies (Marconi et al. 2004; Di Matteo et al. 2005; 
Wang et al. 2006; 2008), there is no doubt that the detection of 
the fuzz will provide the most powerful diagnostic to evolution- 
ary chains of galaxies and quasars from details of the radiative 
feedback. 

Finally, it should be mentioned that X-ray emission from fuzz 
is extremely faint unless the maximum Lorentz factor is low- 
ered by a factor of 10. This could be plausible if the factor A 
decreases by a factor of 10. In such a context, the fuzz is able 
to produce hard X-rays detected by the Nustar (The Nuclear 
Spectroscopic Telescope Array). 

5. CONCLUSIONS AND DISCUSSIONS 

In this Letter, we show that the enhanced radiation pres- 
sure inevitably leads to strong shocks in the dusty gas of both 
homogeneous and clumpy ISM. The first Fermi shock accel- 



eration results in production of relativistic electrons in the 
shocked ISM, from which produces non-thermal emission via 
synchrotron radiation and inverse Compton scattering, display- 
ing multiwavelength fuzz from radio, infrared and 7-ray bands. 
The current radio telescopes, MIPS, HESS and VERITAS are 
capable of detecting the fuzz. The multiwavelength image of 
the fuzz provides an interesting and new clue to understand the 
interaction between quasars and their environments. 

The predicted 7-ray emission from fuzz could have signif- 
icant contribution to cosmic 7-ray background if all quasars 
have to undergo such a red phase during their cosmic evolution. 
We will issue this problem in the future. 
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